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Differential conductance spectra are obtained from nanoscale junctions on the heavy-fermion 
superconductor CeCoIns along three major crystallographic orientations. Consistency and repro- 
ducibility of characteristic features among the junctions ensure their spectroscopic nature. All 
junctions show a similar conductance asymmetry and Andreev reflect ion- like conductance with re- 
duced signal (~ 10% — 13%), both commonly observed in heavy-fermion superconductor junctions. 
Analysis using the extended Blonder- Tinkham-Klapwijk model indicates that our data provide the 
first spectroscopic evidence for d^-2_y2 symmetry. To quantify our conductance spectra, we propose 
a model by considering the general phenomenology in heavy fermions, the two-fluid behavior, and 
an energy-dependent density of states. Our model fits to the experimental data remarkably well 
and should invigorate further investigations. 

PACS numbers: 74.50.-|-r, 74.20.Rp, 74.45. 4-c, 74.70. Tx 



The 1-1-5 family of heavy fermion (HF) compounds, 
CeTIns (T — Co, Rh, Ir), has attracted great interest 
because of their novel and rich physical phenomena [l[. 
Various thermodynamic and transport measurements in- 
dicate that the superconducting order parameter (OP) 
in CeCoIns is d-wave [ij, but precise locations of the 
line nodes over the Fermi surface remain controversial 
0, Q . Not only is the interpretation of these experiments 
complex but also they intrinsically cannot provide 
phase information of the OP. In this Letter, we report 
differential conductance data on CeCoIns as a function of 
temperature, magnetic field, and crystallographic orien- 
tation. Our results show the first spectroscopic evidence 
for dx2_y2 symmetry. We further present a model, which, 
for the first time, quantifies the Andreev signal and con- 
ductance asymmetry observed in normal-metal/heavy- 
fermion superconductor (N/HFS) junctions [1, SSQ- 

Andreev reflection (AR), the scattering of a quasipar- 
ticle off an attractive pair potential, occurs at an N/S 
interface as the retro-reflection of an electron as a hole 
0. If the N and S are in good electrical contact and 
their Fermi velocities are well matched, the conductance 
is twice the normal state value within the superconduct- 
ing energy gap, A, and rapidly returns to the normal 
state value outside A, directly providing energy gap in- 
formation. The Blonder-Tinkham-Klapwijk (BTK) the- 
ory [lol| describes the conductance data in practical N/S 
junctions remarkably well using a dimensionless barrier 
strength parameter, Z^s — [Zq + (1 — r)^/4r]-^/^, where 
Zq is due to a physical potential barrier and r is the ratio 
of Fermi velocities As Z^s increases from zero, the 
junction moves from the AR to the tunneling regime. In 
HF systems, the electronic mass is highly enhanced (by 
~ 10 — 10^) with a correspondingly reduced Fermi veloc- 
ity. According to the above formula, a N/HFS junction is 



inherently in the tunneling regime and AR cannot occur. 
However, AR is frequently observed in N/HFS junctions 
(Ref. and references therein). Deutscher and Nozieres 
addressed this discrepancy by assuming the boundary 
conditions are not affected by mass enhancement jTl] |. 

Three crystallographic surfaces of high quality 
CeCoIns single crystals are prepared: the (001) face of 
as-grown crystals, and the (110) and (100) surfaces by 
embedding into epoxy and polishing. X-ray diffraction 
analysis confirms their crystallographic orientations [§]. 
A light HCl etch exposes fresh surfaces prior to junction 
formation. The average surface roughness, which ranges 
1 — 2 nm over ~ 1 /im^ area, is much smaller than the 
junction size, estimated below. Electrochemically pol- 
ished metal tips, Au (primarily), Al, or Nb, are used 
as counter-electrodes. Nanoscale junctions are formed 
by electromechanical adjustment of the tip-sarnple dis- 
tance in our Cantilever- Andreev- Tunneling rig 12] that 
is run in a ^He fridge. Differential conductance (dl/dV) 
data are taken using a standard four probe lock-in tech- 
nique over wide ranges of temperature (400 mK — 60 
K) and magnetic field (0 — 9 Tesla). High-bias junc- 
tion resistances typically range 1 — 5 f2, which corre- 
spond to junction sizes of 20 — 50 nm, estimated us- 
ing Wexler's formula [l3|, indicating that junctions are 
in the ballistic (Sharvin) limit; the extreme cleanness of 
this compound makes it readily accessible (the electronic 
mean free path ranges several micrometers at low tem- 
perature [lj|). However, we stress this is not a sufRcient 
condition for AR spectroscopy (see Ref. [3] and refer- 
ences therein). Reproducibility and consistency of the 
conductance spectra along different crystallographic ori- 
entations are crucial to ensure their intrinsic and spec- 
troscopic nature. 

Normalized conductance spectra for the (001) and 
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FIG. 1: (color online). Normalized conductance spectra of 
CeCoIns/Au junctions (a) along (001) (after Ref. Q|) and 
(b) along (110) orientations. Data are shifted vertically for 
clarity. Note the temperature evolution of the background 
conductance, whose asymmetry is quantified in (c) by the 
ratio between conductance values at -2 mV and at +2 mV 
in (a); the inset is a semi- logarithmic plot (T* is the HF 
coherence temperature). Junctions along three orientations 
are compared in (d) at ~ 400 mK and in (e) at ~ 1.5 K. 



(110) junctions are displayed in Figs. 1(a) and 1(b), re- 
spectively. At high temperatures, the conductance curves 
of the (001) junction are symmetric and flat; character- 
istic of simple metallic junctions. As the temperature is 
reduced, they become asymmetric and curved. This con- 
ductance asymmetry begins at the HF coherence tem- 
perature T* (~ 45 K) [16| and increases with decreasing 
temperature down to Tc (2.3 K), below which it remains 
constant. The same behavior is observed in the (110) 
junction and the data near and below Tc are shown in 
Fig. 1(b). A plot of the ratio of the conductance values at 
-2 mV and +2 mV quantifies this asymmetry (Fig. 1(c)). 
According to the two-fluid model proposed by Nakatsuji, 
Pines and Fisk [l^ , the spectral weight for the emerging 
HF liquid grows below T* and saturates below Tc |lq . 
and our conductance data track this behavior. 

The conductance near zero-bias begins to be enhanced 
as Tc is crossed and increases with decreasing temper- 
ature, indicating its origin is AR. Conductance data at 
two temperatures are compared for all three orientations 
in Figs. 1(d) and 1(e). Three consistent and reproducible 



FIG. 2: (color online). Comparison of conductance data: (a) 
the (100) and (b) the (110) junctions, (c) and (d) Calculated 
conductance curves using the d-wave BTK model (F = and 
T = 0) for antinodal and nodal junctions, respectively, (e) 
and (f) Magnetic field dependence for the (100) junction at 
400 mK and the (110) junction at 420 mK, respectively. 



characteristics are observed at low temperature, indicat- 
ing we are sampling intrinsic spectroscopic properties. 
First, all spectra are asymmetric with the positive-bias 
side (electrons flowing into CeCoIns) always lower than 
the negative-bias branch. We have seen similar conduc- 
tance asymmetry from more than two hundred junctions 
on pure and Cd-doped CeCoIns along all three direc- 
tions. This is in strong contrast with the symmetric con- 
ductance data we obtained from junctions on non-HFS 
such as Nb, MgB2 [11], and LuNi2B2C. All these obser- 



vations strongly indicate that the conductance asymme- 
try arises from intrinsic properties in CeColns. Second, 
the conductance enhancement occurs over similar voltage 
ranges, ~ ±(1 — 1.5) mV. Third, the normalized zero-bias 
conductance (ZBC) ranges 1.10 — 1.13, showing that our 
observed Andreev signal is much smaller than the theo- 
retical prediction of 100% [l3|. We reported Q that it 
is too small to fully account for the conductance spec- 
tra using the existing BTK models even considering the 
mismatch in Fermi surface parameters, nonzero Zeg, and 
large quasiparticle lifetime broadening factor (F). Our 
model proposed below enables us to quantify it success- 
fully and elucidates properties of the HFS state. 

Conductance spectra for in-plane junctions are plotted 
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in Figs. 2(a) and 2(b). While both spectra exhibit similar 
background asymmetry, differences in the sub-gap region 
are striking. The (100) data appear rather flat, similar to 
the (001) junction, whereas the (110) data are cusp-like. 
This shape difference persists even to higher temperature 
despite an enhanced thermal population effect, indicat- 
ing it is intrinsic. We compare these data with calcu- 
lated conductance curves using the d-wave BTK model 
[13], as shown in Figs. 2(c) and 2(d) for antinodal and 
nodal junctions, respectively. Both curves are identical at 
Zeff = but quickly evolve in dramatically different man- 
ners with increasing Z^s- For an antinodal junction, the 
ZBC is gradually suppressed and a double-peak structure 
develops for Z^s ~ 0.3. For a nodal junction, the ZBC in- 
creases and the sub-gap conductance narrows into a sharp 
peak. This is the signature of Andreev bound states 
(ABS) which arise directly from the sign change of the 



OP around the Fermi surface [18[ . We stress the flat con- 
ductance shape observed in the (100) junction can only 
occur for an antinodal junction with Zes — 0.25—0.30 but 
cannot occur in a nodal junction at any Zee value. Mean- 
while, the cusp-like feature in the (110) junction cannot 
occur in an antinodal junction unless Zcff is small enough 
0.1), an unlikely condition in N/HFS junctions; it 
can only be explained by a sign change of the OP, rul- 
ing out anisotropic s-wave. We therefore assign the (100) 
and (110) orientations as the antinodal and nodal direc- 
tions, respectively, providing evidences for dj,2_j^2-wave 
symmetry and resolving the controversy on the locations 
of the line nodes 0, H 0]- Note this is a spectroscopic 
measurement of the superconducting OP symmetry in 
that it can detect its sign change, in contrast with other 
measurements that probe only the gap anisotropy, in- 
cluding heat transport and NMR. The ABS-originated 
ZBC peaks are reported to split spontaneously and/or 
under applied magnetic field in high-Tc cuprate tunnel 
junctions [l^. One of the widely adopted explanations 
is the Doppler shift of ABS. We test it by applying a 
magnetic field perpendicular to the ab plane, a configu- 
ration for a maximal shift, if any. As shown in Figs. 2(e) 
and 2(f), no splitting but only gradual suppression of the 
Andreev signal is observed in both junctions. We have 
observed that the Z^s values do not change significantly 
by using different tips (e.g., see Fig. 4(b)), samples, and 
contact pressures. While further studies are necessary, 
diminution of the Doppler effect due to large junction 
transparency, small tunneling cone, and atomic-scale dis- 
order, has been suggested to explain similar behaviors in 
some cuprate junctions ^2^. 

Both asymmetric conductance and a reduced An- 
dreev signal have been commonly reported by others in 
CeCoIns fi\ and other HFS junctions The temper- 
ature dependence of the conductance asymmetry (Fig. 
1(c)) qualitatively follows that of the spectral weight of 
the coherent HF liquid in the two-fluid model 16] pro- 
posed to be general to HFs [2l|. These observations 
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FIG. 3: (color online), (a) Normal state conductance spectra 
of the (001) junction, whose features are reproduced qualita- 
tively by the simulated Lorentzian DOS curves in (b). 



strongly suggest a universal mechanism to explain the 
charge transport at HF interfaces. We attribute it to the 
emergent HF liquid: with decreasing temperature, the 
logarithmic increase of the electronic specific heat coeffi- 
cient signals the increase of the electronic mass and, 
cquivalcntly, the electronic density of states (DOS) at 
the Fermi level. The electrical conductance of a clean 
metallic junction, where one electrode is a simple metal, 
is given by: |f (F) cx JJ vD{e)^^^j^^dedn, where v 
is the velocity, D(e) the DOS of the counter-electrode, 
/ the Fermi-Dirac distribution function, and dQ the dif- 
ferential solid angle Q. For a simple metal, the DOS 
is constant around the Fermi level and thus divides out 
in the normalized conductance giving flat and symmetric 
shape. For a HF metal, the energy-dependent DOS is re- 
flected in the conductance data [221; an asymmetric DOS 
yields asymmetric conductance. In Fig. 3(a), our normal 
state conductance data exhibit an increasing asymmetry 
with decreasing temperature while the conductance peak 
sharpens and shifts towards the Fermi level. We have 
investigated several proposed models in reproducing this 
non-trivial temperature dependence, including large See- 
beck effect and non- Fermi liquid behavior in HFs. We 
find that only our model, in which we assume a peak in 
the DOS below the Fermi level, reproduces the observed 
experimental features (Fig. 3(b)) |23| . 

Measurements of de Haas-van Alphen effect in 
CeCoIns show that heavy and light fermions reside on 
disparate Fermi surfaces 



24l |. Tanatar et al. reported 



[25l | that the light electrons remain uncondensed below 
Tc. The junction conductance would therefore be a mea- 
sure of both normal and superconducting channels. We 
modify the BTK model such that the total conductance 
is given by the sum of two parallel conductance chan- 
nels: ^{V) ^ cu^,jL\^{V) + (1 ~ u;^,)^\m, where the 
first term is the usual BTK conductance [l7| due to the 
superconducting heavy electrons, the second a constant 
conductance due to the uncondensed light electrons, and 
ujh the weighting factor related to the HF spectral weight 
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FIG. 4: (color online), (a) Best fit to the data using our 
modified BTK model, (b) Comparison of the data with cal- 
culated DOS curves. Right inset: Fano lines simulated by 
/(e) — (qp + e)^/(l + e^), where qf is the Fano factor and 
e normalized energy. Left inset: representative conductance 
data of Fano line shape taken with Al tips. 



7] and A are the peak height and half width, respectively. 

Figure 4(a) shows the best-fit curve for the (001) data 
taken at 400 niK with w,, = 0.51, A = 600 ^eV, F = 95 
/xeV, Zcff = 0.28, ?7 = 1, A = 5 meV and eo = -2.1 meV. 
The physical origin of the peak in the DOS below the 
Fermi level remains an open question, particularly with 
direct measurements such as photoemission or tunneling 
still lacking. Compared to our previous fit @ using the 
single channel BTK model 17|, the quality of the fit is 
remarkable; it nicely reproduces both features of a re- 
duced AR signal and the conductance asymmetry. The 
obtained energy gap gives = 6.05, suggesting 

strong coupling in agreement with literature. We find the 
fit using our proposed model is particularly sensitive to 
LOh', for smaller values it becomes much poorer; for larger 
values it requires smaller A and larger F values, causing 
the same problem of an unphysical temperature depen- 
dence of F as in our previous analysis [6[. Thus, both 
superconducting and normal conductance channels are 
necessary. Our model provides a natural explanation for 
the unreduced Andreev signal in our CeCoIns/Nb junc- 
tions at temperatures between the two T^'s [2^; here 
both heavy and light electrons in CeCoIns participate 
in AR. In the high temperature region, the fitted curve 
deviates substantially from the data and we find that a 
Fano 27| line shape as a background gives a much bet- 
ter fit, as illustrated in Fig. 4(b). The Fano line as a 
function of the Fano factor, gp, is shown in the right in- 
set of Fig. 4(b). The line shape for q-p ^ -2 is similar 
to those observed in various junctions by GoU et al. 
and by us over wide voltage ranges, as exemplified in the 
left inset of Fig. 4(b). Since the physical origin of the 
Fano-like background needs to be clarified further, we 
defer quantitative analyses of the complete conductance 



spectra [23|. As seen, our model correctly captures the 
underlying physical phenomena. This implies that micro- 
scopic details such as boundary conditions do not play a 
key role in reducing an Andreev signal. Since the two- 
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fluid behavior is proposed to be universal to HFs 
our model should also be generally applicable. 
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